Work extraction in the spin-boson model. 
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We show that work can be extracted from a two-level system (spin) coupled to a bosonic thermal 
bath. This is possible due to different initial temperatures of the spin and the bath, both positive 
(no spin population inversion) and is realized by means of a suitable sequence of sharp pulses applied 
to the spin. The extracted work can be of the order of the response energy of the bath, therefore 
much larger than the energy of the spin. Moreover, the efficiency of extraction can be very close 
to its maximum, given by the Carnot bound, at the same time the overall amount of the extracted 
work is maximal. Therefore, we get a finite power at efficiency close to the Carnot bound. 

The effect comes from the backreaction of the spin on the bath, and it survives for a strongly 
disordered (inhomogeneously broadened) ensemble of spins. It is connected with generation of 
coherences during the work-extraction process, and we derived it in an exactly solvable model. All 
the necessary general thermodynamical relations are derived from the first principles of quantum 
mechanics and connections are made with processes of lasing without inversion and with quantum 
heat engines. 

PACS numbers: PACS: 03.65.Ta, 03.65.Yz, 05.30 



1. INTRODUCTION. 



A known feature of technological progress is the increase of human ability to control and design the microscopic 
world. Recent efforts in manipulating simple quantum systems e.g., in the context of quantum computing or quantum 
chemistry, is one aspect of this general trend. Another aspect is the field of quantum thermodynamics whose main 
objective is in designing and studying new thermodynamical processes in the domain where quantum features of 
matter are relevant. In particular, this activity aims to improve our understanding of the standard thermodynamics 
P, 0, Q by addressing its concepts from the first principles of quantum mechanics 0, 0, ■ The current activity in 
quantum thermodynamics includes quantum engines ja El El E3 , general aspects of work extraction from 
quantum systems |l4J , thermodynamical aspects of quantum information theory |15l [lfj , and limits of thermodynam- 
ical concepts such as the second law 0, 0, 0] and the temperature fl?|. There were also much earlier applications 
concerning, in particular, thermodynamic aspects of lasers and masers [l8| . 

Our prese nt purp ose is to study work extraction from a two-temperature system on the basis of the known spin- 
boson model |20l l2ll l22l l23| : a two- level system coupled to a bosonic thermal bath. The motivation to use a two-level 
system is nearly obvious, it is almost everywhere, and it is the minimal model having non-trivial quantum features. 
The necessity of the bath has to be stressed separately since, in the usual practice of quantum systems manipulation, 
the bath is a serious hindrance. As follows, the process of work-extraction really needs external thermal baths: The 
second law in Thomson's formulation — which is derived as a theorem in quantum mechanics |26l 1271 1281 l29i l3Ct l31| — 
forbids work-extraction from an equilibrium system by means of cyclic processes generated by external fields. The 
easiest way to employ an equilibrium system in work-extraction is to attach it to a thermal bath having a different 
temperature, thus forming a local-equilibrium state. The overall system is then out of equilibrium and work-extraction 
from cycles is not forbidden, at least in principle. This was shown explicitly in 0,0,1a- 

This general restriction determined the way how standard quantum work-extraction (also known as amplification 
or lasing/masing) processes are designed [l8|. The most traditional lasers and masers operate by extracting work 
from an ensemble of two-level systems having a negative temperature, in other words, population inversion, which is 
a strongly non-equlibrium state. More recent schemes of lasing without inversion employ non-equilibrium states of 
three (four, multi) level systems without population inversion of energy levels, but with initially sizable non-diagonal 
terms of the corresponding density matrix in the energy representation, usually called coherences |32l l33l |. These 
schemes attracted attention due to both their conceptual novelty and the fact that non-zero non-diagononal elements 
represent a weaker form of non-equilibrium than population inversion, and thus their preparation can be an easier 
task |3l|33. 

The mechanism of work-extraction proposed in the present paper differs from the standard ones in several aspects: 

• Work-extraction (amplification, lasing) can be achieved in two-level systems without population difference and 
without using an initially coherent state. A setup consisting of a positive temperature spin interacting with a 
thermal bath at some higher or lower temperature suffices to extract work and thereby amplify pulsed fields 
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acting on the spin. Moreover, the extracted work can be of the order of the bath's response energy, which is 
larger than the energy of the spin. Thus, when viewed as lasing without inversion, the presented mechanism 
offers definite advantages over the existing schemes. 

• The effect survives for a disordered ensemble of spins, where the spin have a random energy with a large 
dispersion. The reason of the survival is the possibility to combine the work-extraction process with the spin-echo 
phenomena |34Ll35j]. As a consequence, we have a phenomenon even more amazing than the original spin-echo: 
a high-temperature, completely disordered ensemble of spins can serve as a medium of work-extraction. 

• The efficiency of work-extraction can approach its maximally possible value given by Carnot bound. Moreover, 
the efficiency is maximized simultaneously with the overall amount of the extracted work. In addition, the power 
of work (i.e. the work divided over the total duration of the work-extraction process) is finite. Thus, in marked 
contrast to the original Carnot cycle 0, |j| and some of its realizations in quantum engines 0] 1 the three basic 
objectives of a good work-extraction process are met together: large amount of extracted work, high efficiency, 
and finite power. 

The origin of the presented mechanism is that, besides well-known effects of dissipation and decoherence induced 
by a thermal bath on a spin interacting with it, there is another effect the presence of which is frequently not 
acknowledged. This is the backreaction of the spin to the bath, which in combination with external fields influences 
the spin's dynamics. The effect exists even for relatively small — but generic — bath-spin coupling constants, but is 
typically neglected from standard weak-coupling theories [l9j |. Our present treatment of the bath-spin interaction is 
exact and allows to study the full influence of backreaction and memory effects. 

This paper is organized as follows. In section[21we recall a version of the spin-boson model we work with. It nowdays 
became one of the most popular models in the theory of open quantum systems [2(3, COl H3, H3, H3 ■ Section[5]describes 
the action of external fields and discusses the definition of work. The next section presents experimental realizations 
of our setup in various situations of two-level system(s) interacting with a thermal bath. Section [5] discusses general 
limits of work-extraction from a two-temperature system. The next two sections describe our basic results on the 
work-extraction, effciency, and the power of work. The last section offers our main conclusions and compares our 
results with the ones existing in literature. Several technical questions are considered in Appendices. We have tried 
to make this paper reasonably self-contained. This especially concerns the concepts and relations of the standard 
thermodynamics, which are not accepted uncritically, but, in many situations, are derived from the first principles of 
quantum mechanics. 



2. THE MODEL 

As common when dealing with open systems, the Hamiltonian H is composed by three parts: 

H = H a +H B + H 1 . (2.1) 
Hs stands for the Hamiltonian of a two level system (spin |): 

Hs = ~&z, £ = m, (2.2) 

where a x , a y and a z are Pauli's matrices, and where the energy levels are ±|. 

The spin interacts with a thermal bath which is a set of harmonic oscillators. In some cases this may be taken in 
the literal sense, when harmonic oscillators represent phonons or photons. It is also known that rather general classes 
of thermal baths can be effectively represented via harmonic oscillators 0, 0] ■ Thus for the Hamiltonian of the 
bath we take: 

Hb = ^ fuv k a\a k , [au a' k ] = Su, (2.3) 

k 

where at and &k are creation and annihilation operators of the bath oscillator with the index k. The thermodynamic 
limit for the bath will be taken later on. 

The next important point is to specify the interaction between the spin and the bath. Recall that any reasonable 
model of a thermal bath is expected to drive a non-stationary state of the spin towards a stationary state. In this 
res pec t, for two-level systems, one distinguishes two types of relaxation processes and the corresponding times scales 

i III EE El & 
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1. 72-time scale related to the relaxation of the average transversal components (a x ) and (<r y ) of the spin (de- 
coherence). Note that the very notion of the transversal components is defined by the form (|2.2|) of the spin 
Hamiltonian. 

2. 7i-time scale related to the relaxation of (<t z ). It is customary to have situations, where 

T 2 «Ti, (2.4) 

the main physical reason being that the transversal components are not directly related to the energy of the 
spin. 

Our basic assumption on the relaxation times is 1|2.4[) 1 . Moreover, to facilitate the solution of the model we will 
disregard Ti time as being very large, thereby restricting the times of our interest to those much shorter than T\. The 
interaction Hamiltonian is thus chosen such that it induces only transversal relaxation: 

Hi = ^Xa z , X = J29k(a{ + a k ), (2.5) 

k 

where g k are the coupling constants to be specified later, and where X is the collective coordinate operator of the 
bath. 

The last ingredient of our model is external fields which are acting on the spin. However, before discussing them in 
the next section, we shall recall how the model with Hamiltonian Eq. (|2.1I) is solved without external fields. 

1. Heisenberg equations and their exact solution. 

Hcisenberg equations for operators a z (t) and a k (t) read from l|2.2l 12.31 12.51 12. 1|) : 

& z =0, a 2 (t) = a z (0), (2.6) 



Eqs. (E21 E2J are solved 



and then 



where 



a k = -[H,a k \ = -iuj k a k - -g k a z . (2.7) 



a k {t) = e-^a k (0) + (e~^* - l) , (2.8) 

2u k 



X{t)=f,{t)-& z G{t), (2.9) 



Git) = V^-(l-coso; fc i), (2.10) 

u ^ k 
k 

quantifies the reaction of the spin on the collective operator of the bath, and where we denoted 

fi{t) =$>*[4(0K Wfe < + a fe (0)e-^ 4 ], (2.11) 
fc 

for the quantum noise operator 2 . Recalling the standard relations 

<7± = a x ± ia y , [6- z ,a±] = ±2a±, a z a± = ±a ± , (2-12) 



1 There is also a third relaxation time which has a different origin. It only appears when dealing with an ensemble of non-interacting 
spins each having Hamiltonian 12.21 with a randomly distributed energy e (dephasing). The influence of is studied in section 171 

2 Note that the commutator of the quantum noise is a c-number: [r)(t), fi(s)] = —2i ^2 k sinLjj,(t — s) = —2ihG(t — s). 



4 



and using i|2.9[l and [X(t), a±(t)] = — since they belong to different Hilbert spaces, — one derives 

fr± = j[H,(t±] = ±i (o + i") a+=i (±Cl±fj(t) - G(t)) a ± . 
These equations are solved as: 

a±(t) = exp[±mt-iF(t)] n±(0, t) <t±(0), 



U±(t ,ti) = T exp 



±i 



ds fj(s) 



F(t) 



sin uj^-t 



(2-13) 

(2-14) 
(2-15) 

(2-16) 



where T stands for the time-ordering operator. It is seen again from l|2.I4[l that there are two effects generated by 
the bath-spin interaction: besides random influences entering with the quantum noise T}(t), there is a deterministic 
influence generated by the backreaction term F(t), somewhat similar to damping (friction) in the problem of quantum 
brownian motion. 



2. Factorized initial conditions. 



Let us assume that initially, at the moment t = 0, the bath and the spin are in the following factorized state: 

e -/3H B 

p(O) = p s (O)®p B (O) = ps(O)0 -jp (2.17) 

where ps(fi) is the initial density matrix of the spin, and where the bath is initially at equilibrium with temperature 
T=l/f3. 

Factorized initial conditions are adequate when the spin is prepared independently from the equilibrium bath and 
then is brought in contact to it at the initial time 3 . For example, injection of an electronic spin into a quantum 
dot, or creation of an exciton by external radiation. Yet another situation where factorized initial conditions can be 
adopted is a (strong) selective measurement of a z by an external apparatus. In this case ps(0) is an eigenstate of a z 
upon which the selection was done. Non- factorized initial states are commented upon below, in section \2. 41 

The equilibrium relation 

(4(0)) = <a*(0)> = 0, (2.18) 
(4(0)a*(Q) + a fc (0)4(0)) = coth (2.19) 

derived from l|2.17|l . imply that the quantum noise is a stationary Gaussian operator with 

m)) = 0, (2-20) 

and having the time-ordered correlation function: 

K T (t-t') = (T[fj(t)fi(t')})fi =J2 9k coth^^p^ cosu k (t-t')- isgn (*-*') sin u> k (t-t') (2.21) 

where the average (...) is taken over the initial state (|2.17l) . It can be written as 

K r {t) = K{t)-iG{t), (2.22) 



3 It is useful to note that this process of bringing spin in contact to the bath need by itself not be connected with any fundamental energy 
cost. Imagine, fo example, a sudden switching of the interaction Hamiltonian Hi = ^Xa z . Since in the equilibrium state of the bath 
(X) = 0, the work done for the realization of this switching is zero. 
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where 



K(t - f ) = KK T (t - t') = \{mv{t') + ^0<K*)) = E 9l coth cos Wfc (f, - f), (2.23) 

is the symmetrized correlation function. 

Since ?)(t) is a gaussian random operator, one can use Wick's theorem for decomposing higher-order products 4 . 
Due to the factorized structure (|2.17() of the initial state, the common averages of rj and various spin operators can 
be taken independently. For example, averaging Eq. I|2.14[l and using Wick's theorem together with the arithmetic 
relation k\ 2 k (2k - 1)!! = (2fc)! one gets: 

(M*)) = e ±mt - lF{t) ( n±(0, t) } (<t ± (0)) = e ±int ~^ <<t±(0)), (2.24) 

where for t% > t\: 



\ °° f— \\ k Z"* 2 Z"* 2 
^/^W / "7 d Sl ...d S2fe (T[rK S i)...r)( S2fc )]) 



fc=0 

= exp 

and where 



ti "'ti 

= exp[-£(t 2 -t 1 )+iF(i 2 -t 1 ) 



2 J J ds i ds 2^r(si - s 2 ) 



(2.25) 



£(*) = - / / d Sl d S2 K( Sl - S2 ) = f dsi f 1 ds 2 K(s 2 ). (2.26) 
^ Jo Jo Jo Jo 

As seen from (I2.24|) . characterizes the decay of (a±) due to the interaction with the bath. 

3. Ohmic spectrum of the bath. 

The coupling with the bath can be parametrized via the spectral density function J(ui): 

J{u)=Y J 9l5{"-"k). (2.27) 

k 

In the thermodynamical limit the number of bath oscillators goes to infinity, and J(oj) becomes a smooth function, 
whose form is determined by the underlying physics of the system-bath interaction. 
We shall be mainly working with the ohmic spectrum: 

J(w)=7wr u/r . (2.28) 

where 7 is the dimensionless coupling constant, and where T is the maximal characteristic frequency of the bath's 
response. This spectrum and its relevance for describing quantum open systems was numerously discussed in literature; 
see, e.g. 



4 Recalling Wick's theorem (for (77) = 0): Any correlation of an odd number of »7's vanishes. A correlation of an even number of fj's is 
equal to the sum of products of pair correlations, the sum being taken over all pairings. For example: 

(T fj(ti)v(t2)v{t 3 )v(U)) = {T f)(ti)v(t2))(T fj(t 3 )fj(U)} + <T ft(ti)fi(t 3 )}{T fi(t 2 )fi(t 4 )) + (T fj{ti)f,(t 4 )){T f,(t 2 )f,(t 3 )) . 

Note that the similar Wick-decomposition of (T f)(ti) . . .»j(t2fc)} will be a sum of (2k — 1)!! = (2k — l)(2k — 3). ..3 terms. Wick's theorem 
is related to the fact that the commutator of the quantum noise is a c-number; see Footnote 121 
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a. Quantum noise correlation function and decay times. 



The correlation function of the quantum noise in the ohmic case, using is given by : 



K(t) 



dto J{ui) coth 



2T 



/>oo 






/ divu! coth 


e 


10 


2T_ 





cos cot 
~" /r cosut. 



(2.29) 
(2.30) 



Recall that the decay factor is related to K(t) via Eq. I|2.2(j[) : £(t) — K(t). Properties of these functions are 
worked out in Appendix 1X1 In particular, for one gets from Eqs. (|A13I lA14l) the following exact expression: 



^(t) = 7 In 



r 2 (i + _r) yiTW 
r(i + ^-^)r(i + ^ + if 



(2.31) 



where T is Euler's gamma function. It is seen that the temperature is controlled by the dimensionless parameter 
T/(hT). 

Let us now determine the behavior of this quantity for low and large temperatures. Using Eq. (|A7() one obtains for 
hT/T ^> 1 (low temperatures): 



m = 7 In 



h(3 , firt 
— sinn — — 
irt \h(3 



This implies two regimes of decay: power-law and exponential. 

t«fi/3: e-«W = (l + r 2 t 2 )"^ 2 , 



-f(t) 



For HT/T <C 1 (high temperatures) one uses Ea. (|A12|l to get 



2 7 T 

"fir" 



rtaxctan(n) - - ln(l + T 2 t 2 ) 



(2.32) 

(2.33) 
(2.34) 

(2.35) 



This time the possible regimes of decay can be approximated as gaussian and exponential. 



t> i/r 



/j 2 



7 Tr' 



e -f(t) = e -t/r a % 



2jT 



(2 


36) 


(2 


37) 


with 


the 



In this latter case, as seen from Eq. (|A12(I . K(t) behaves as approximate delta-function: K(t) ~ n(i+t^r^) 
strength 2jT/h determined by parameters 7 and T. Note that in all the above cases the characteristic times of decay 
become shorter upon increasing the temperature T or coupling constant 7, as is expected. The gaussian regime of 
decay was also numeruosly observed in NMR experiments (see [44J and refs. therein). This regime is the basis of the 
quantum Zeno effect [2(j and was recently predicted to govern the reduction process in quantum measurements . 



b. The G-factor. 

Finally we will indicate the form of the backreaction functions G(t) and F(t) in the ohmic case (see Ea. (|2.10|l 'l. As 
will be seen below, these functions are rather important for our purposes. 



1 



F(t) = 7 [Tt - arctan(rf)] 



(2.38) 
(2.39) 



Since G{t) becomes equal to a constant on the characteristic time 1/r, it is justified to call the latter the response 
time of the bath. 
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4. Correlated initial conditions. 



Most papers on the system bath models assume factorized initial conditions. However, in many situations the use 
of such a condition is difficult to justify a priori, since it implies a possibility of switching the system-bath interaction. 
Non-factorized initial conditions can be Gibbsians that are modified at the initial time, as considered in 6] for the 
Caldeira-Leggett model and in 5] for the spin-boson model. 

For our present purposes it is sensible to use the following correlated initial conditions for the spin and the bath: 



P(0) 



■ exp 



-I3 S H S ~ /3(Hi + H B ) 



Z = tr e 



-/3 s B s -/3(Hj+H B ) 



(2.40) 



where f3 is the inverse temperature of the bath and /3s is that of the spin. 

The initial condition (|2.40l) with (3g ^ f3 can be generated from the equilibrium equal-temperature state of the 
overall system via cooling or heating the bath by means of some superbath. During this process a z is conserved, and 
the bath relaxes to its new temperature under an "external field" ± \X generated by the interaction Hamiltonian Hi 
with a z = ±1. More details of this procedure are given in AppendixIBl 

In the thermodynamical limit for the bath, the correlated initial condition l|2.40[l is equivalent to the factorized 
condition 1|2.17(1 with 



ps(0) = 



1 



tr e-P$ Hs 



-f3 S H s 



(2.41) 



that is, when starting from the factorized initial condition (|2.17l l2~4l"J) . the dynamics of the overall system builds up 
a correlated state which at times t much longer than the response time of the bath, i > 1/r 5 (ergodic limit), is 
equivalent to i|2.40|l . By saying "equivalent" we mean that the initial conditions i|2.17l I2.41fl and (|2.40() produce the 
same values for spin's observables and for collective observables of the bath (i.e., the ones involving summation over 
all bath oscillators). This equivalence is further discussed in Appendix [HI 

As for the initial state of the spin, it can be deduced from Eq. (|2.40() or from Eq. 12.41(1 



= — tanh 



= 0. 



(2.42) 



In the following we will use the factorized initial condition l|2.17|l since it is technically simpler. The time limit 
t — ► oo will be taken before any perturbation acts on the system to ensure the equivalence with the correlated initial 
condition (|2.40l) . 



3. PULSED DYNAMICS. 



The external fields acting on the spin are described by a time-dependent Hamiltonian 

M*) = \ E h k (t)& k , (3.1) 

k— x,y ,z 

with magnitudes /ifc(i), which is to be added to H defined in (|2.1|l such that the overall Hamiltonian is time-dependent: 

H(t) = H + H F (t). (3.2) 

Eq. Il.'i. l[l represents the most general external field acting on the spin. We shall concentrate on pulsed regime of 
external fields which is well known in NMR and ESR physics [l8l l34l l35l l36l l37l l38l |3t|. For example, it was used 
to describe spin-echo phenomena [34L l35| or processes that switch off undesired interactions, such as those causing 
decoherence pi I2I IM E% . 

A pulse of duration S is defined by sudden switching on the external fields at some time t > 0, and then suddenly 
switching them off at time t + 8. It is well-known that during a sudden switching the density matrix does not change 
0, while the Hamiltonian gets a finite change. Let us for the moment keep arbitrary the concrete form of external 



5 In a more general model, the one considering T\, this limit is in fact t X/T and t <C Ti- 
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fields in the interval (t,t + 5). The Schrodinger evolution operator of the spin+bath from time zero till some time 
t + r, t > S, reads: 



Texp 



t + T 



dsH(s) 



-i(t+T-t-5)6 /hq- 



exp 



t+S 



dsH{s) 



-itH/h 



r "/ R U P (t) e- ltfl ' h . 



(3.3) 
(3.4) 



The LHS of Eq. (|3.3|l contains the full time-dependent Schrodinger-representation Hamiltonian H(s), while in the 
RHS of this equation we took into account that the actual timc-dcpcndcncc is present only between t and t + S. The 
terms er ltE -l n and Q-^+^-t-S^H/H stand for the free (unpulsed) evolution in time-intervals (0,i) and (t + S,t + t), 
respectively. In Eq. I|3.4[) we denoted 



U P (t) 



e^Texp 
Texp 



t+s 



dsH(s) 



t+s 



Texp 



ds e is "l h H F {s + t)< 



sH/h 



(3.5) 
, (3.6) 



for the pulse evolution operator. The transition from 1)3. 5fl to (|3.(i[) can be made by recalling that H(t) = H + Hp(t) 
and then by noting that the expressions in these equations satisfy the same first-order differential equation in S with 
the same boundary condition at S = 6 . 

We focus on pulses so short that the influence of the spin Hamiltonian Mla z /2 and the interaction Hamiltonian 
Hi can be neglected during the interval S. This means that one can take the first term in the Taylor-expansion 
(0 < s < 8): 



e ls "/ h H F {s + t) e 



isH/h 



H F {s + t) 

H F (s + t) 
H F (s + t), 



[H,H F (s + t)] + ... 
— a z + Hi,H F (s + t) 



Thus, for the pulse evolution operator one gets 

U P {t) = T exp 



dsH F (s + t) 



(3.7) 



(3.8) 



The generalization of the evolution operator l|3.4l) to an arbitrary number of short pulses is straightforward. 

Note that in obtaining (|3.8|l we do not require that the bath Hamiltonian H-q during the pulse is neglected. 
Since the external fields are acting on the spin only, the influence of the bath Hamiltonian disappears by itself 
from e lsH / n H F (s + t) e - tsH / h ^ arK } is perfectly kept in the general evolution operator (|3.31 13. 4|) . once the interaction 
Hamiltonian H\ has been neglected. 



In full detail: 



e iSH/h Tcxp 



• / ds H(s) =-{H- e iSH ' h H(t + 8) e ~ iSH / h ) e iSH / h Texp 
Jt J h V / 

= _i e «£/» Hp(t + 8) e -«A/R e«ff/» Texp 



n 

t+S 



t+6 



dsH(s) 



,18 



Texp 



Ase isH/h H ¥ (s + t)e- 



stt/h 



JSH/h + e -iSH/h T 



exp 



ds e iail ' h H F (s + t) e - isH ' % 
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Recalling the orders of magnitude Ml and HjT of the spin energy and the interaction energy, respectively, — in 
particular, recall Eq. (|2.9[) . HG — h J °° dui J(uj)/ui and Eq. I|2.28[) . — one gets the following qualitative criteria for the 
validity of the short pulsing regime 

5<$:mm(n-\ [jT]- 1 ). (3.9) 

As it should be, for very small 7 and a fixed T, the second restriction on 5 is weaker than the first one. More 
quantitative conditions for the validity of the pulsed regime were studied recenlty in the context of decoherence 
suppression by external pulses |24j . 

To deal with the pulsed dynamics in the Heisenberg representation, one introduces the following superoperators: 

£ t A = e ikt ' h Ae- iM ' h , (3.10) 

v t A = ul(t)Au P (t) (3.11) 

Then the Heisenberg evolution of an operator A corresponding to Eqs. I|3.41 13. 8f) reads 

A(t + T) = E t Vt£ T A = e l " t/h Ul{t) e l " T/h Ae~ l " T/h U v {t) e~ l " t/h . (3.12) 



1. Definition of work. 



The action of external fields on the system is connected with flow of work. The work done in the time-interval (0, t) 
is standardly defined as the increase of the average overall energy of the spin and bath defined by the time-dependent 
Hamiltonian ff(f) 0,11 Q: 

W{0,t) = tr[p(t)H(t)} -tr[p(0)-ff(0)]. (3.13) 

Due to the conservation of energy of the entire system (spin+bath+work-source) , work is equal to the energy given 
by the corresponding work-source (source of external fields). 

Since the external fields are acting only the spin, there is a differential formula for the work which uses only 
quantities refering to the local state of the spin and which thus illustrates that the work-sources exchange energy only 
through the spin: 

dW ( , .dH F (t)\ ,„ 
^=tr(^)-Aij, (3.14) 

where Hf(£) as defined by l|3.1f) is the contribution of the external fields into the spin's Hamiltonian, and where ps{t) 
is the density matrix of the spin. Eqs. (|3.13l ETT^f) relate with each other by the von Neumann equations of motion 

p = l -[H{t), p(t)} = l j\H + H F (t), p(t)} (3.15) 

for the common density matrix p(t) of the spin and the bath, where H is the Hamiltonian without external fields 7 . 



7 In order to get 13.131 from 13.141 . note that the external fields are acting only on the spin and dtH-p{t) = dtH(t). Then in expression 
Eq. 13.141 . we can change the reduced density matrix pg for the full density matrix p since the only time dependence of the Hamiltonian 
lives in the Hilbert space of the spin. Then Eq. 13.141 can be written as 

_=tr^)^lij=tr^)^iij. 

Now integrate this expression from to r: 

f o dt^- = W(0,r) = £ dttr L(t)^l\ =tr(p(T)£(T)) - tr (p(O)H(O)) - £ dttr (tft) H(t)) . 
Note that the last integral is equal to zero due to the equation of motion 13.151 . 
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More specifically, we are interested in the work due to a pulse. For the above example of a single pulse at time t 
this quantity reads from (|3.4I 13.121 I3.13[l : 

W{0, t + S) = W(t, t + S)= tr[ (p(t + 5)- p{t)) H] = tr(p(t) [P t H — H]). (3.16) 

This expression is directly generalized to several successive pulses: assume that the pulse V t at time t was followed 
by another pulse Vt+ T at time t + r with r > 0. The work done during the first pulse is given by (|3.16|) . while the 
work done during the second pulse reads: 

W(t + t, t + t + 5) = tr[ (p(t + r + S)-p{t + r)) H] = tr(p(t + r)[Pt+ T H - H}) 

= tr(p(0) £ t V t £r [V t+T H - H]). (3.17) 

Summing this up with W(0, t + S) one gets for the complete work for the 2 pulse situation: 

W(0, t + r + 5) = tr(p(0) [£ t V t £rVt+rH - H}), (3.18) 

as should be. 



2. Parametrization of pulses. 

As seen from (|3.1I I3.ll)l . and taking into account the condition (|3.9|) which, to all effects, can be taken as S — > 0, 
any pulse corresponds to the most general unitary operation in the Hilbert space of the spin (this would correspond 
to a rotation in the classical language). It is convenient to parametrize pulses by coefficients c a ^ as: 

Va a = Ul(t)a a U F (t) = c <*,bh, a = ±,z. (3.19) 

b=±,z 

For more detailed applications we will need the explicit form of Up(t) (see (|3.8ll 'l as a 2 x 2 unitary matrix whose 
determinant can be taken to be unity without loss of generality: 

\ ( ^_( e-^cosi? -e-^sint? 
where 

< 0, tp < 2tt, 0<i?<~. (3.21) 

Parametrizations similar to (|3.20(l are frequently applied in NMR and ESR experiments 0, 13 HE HE HE HE H3 
where the spin is rotated certain degrees over a well-defined axis by tuning the parameters of the laser (microwave) 
pulse applied. 



e 



U p(t) - I ^ sin ^ e ^cos^ ' ' (3 - 20) 



4. REALIZATIONS OF THE MODEL. 



Once the model with all its ingredients has been defined, we discuss some of its realizations and provide some 
numbers. A two-level system coupled to a thermal bath is a standard model for practically all fields where quantum 
systems are studied: NMR, ESR, quantum optics, spintronics, Josephson junctions, etc. Two particular conditions 
are, however, necessary to apply the particular model we study: the condition 1\ 3> T~2 on the characteristic relaxation 
times and the availability of sufficiently strong pulses. On the other hand, we can allow for rather short times 7J, 
since as we will see this timescale can be overcome with the spin-echo technique. 

There are experimentally realized examples of two-level systems which have sufficiently long T2 times, satisfy in 
T\ 3> T2, e.g. 7~i exceeds T2 by several orders of magnitude, and admit strong pulses of external fields. For atoms 
in optical traps, where T2 ~ Is, 1/r ~ 10 _8 s, there are efficient methods for creating non-equilibrium initial states 
and for manipulating atoms by external laser pulses 48]. For an electronic spin injected or optically excited in 
a semiconductor, T2 ~ 1 ps 46], and for an exciton created in a quantum dot T2 ~ 10 _9 s [43; in both situations 
1/r ~ 10~ 9 — 10 _13 s, and femtosecond (10 _15 s) laser pulses are available. In the case of NMR physics T 2 ~ 10~ 6 — 10 3 s, 
1/r ~ 1 ps, and the duration of pulses can vary between lps and 1 ^s HEHEE3- 
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In all above examples the response time 1/r of the bath is much shorter than the internal time 1/0 of the spin. 
Sometimes it is argued that such a separation is related to the large size of the bath and is something generic by 
itself. This is clearly incorrect, since as seen from the derivation in sectional the dimensionless parameter Q/T has 
to do with the form of the bath-spin interaction, rather than with the size of the bath. Moreover, several examples 
of bath-spin interaction are known and were analyzed both experimentally and theoretically, where Q/T ~ 1. For 
example, Ref. [43] focusses on relaxation of nuclear spins with hyperfine frequencies f2 ~ 700 MHz, 72 < 90 MHz, and 
the ratio f2/T may vary between 10 and 0.1. 

Another important parameter that characterizes our setup is the initial polarization | (a z ) | of the spin. It is known in 
NMR and ESR physics that the response of magnetic atoms (nucleus) to external dc magnetic field is best characterized 
by the flC g"°^ cy ratio [lq. which is for example equal to 42 MHz/T for a proton. For an electron this ratio is 
10 3 times larger due to the difference between atomic and nuclear Bohr magnetons, and for 15 N it is 10 times 
smaller. Thus at temperature T — IK and magnetic field B = IT the equilibrium polarization of a proton is only 
|((7 Z )| = tanh 2k^r = 10~ 3 7 while for an electron it is ~ 1. 



1. Exact solution versus various approximations. 

The model as stated above that is, with the Hamiltonian I|2.1II3.2|1 — is exactly solvable for all temperatures and 
all bath-spin coupling constants. It is useful at this point to recall the reader what are the specific reasons to insist 
on this feature. The model with Hamiltonian 1|2.1[1 is a particular case of a more general spin-boson model, where the 
influence of 7i-time is retained either via an additional term oc a x in the Hamiltonian of the spin, or via an additional 
coupling in the interaction Hamiltonian. This model is in general not solvable, and what is worse there are no realiable 
approximate methods which apply for a fixed (maybe weak) coupling to the bath and for all temperatures including 
the very low ones. The standard weak coupling theories — both markovian leading to well-known Bloch equations, 
and non-markovian ones — are satisfactory only for sufficiently high temperatures, while at low-temperatures weak- 
coupling series are singular, and different methods of their resummation produce different results. In this context, 
compare, e.g., convolutionless master equations extensively discussed in 20] with a convolutional one worked out in 

m 

This situation becomes even more problematic under driving by external fields. The objects studied by us — such 
as work, energy of the spin — can be rather fragile to various not very well-controlled approximations, since there 
are general limitations governing their behavior: Thomson's formulation of the second law and restrictions on work 
extraction from a two-tem pera ture s y stem (discussed below). These limitations are derived from the first principles 
of quantum mechanics |26l. 1271 12^ . f29t l3tX l31| and have to be respected in any particular model. 



5. GENERAL RESTRICTIONS ON WORK EXTRACTION. 



The setup of two systems having initially different temperatures and interacting with a source of work allows to 
draw a number of general relations on work-extraction. Starting from the following general assumptions: 

1. Out of equilibrium initial conditions. The initial conditions at the moment t = are given by Eq. (|2.40|) . where 
the bath and the spin have initially different temperatures T and Ts, respectively. Recall from discussion in 
section |2~H that after a small lapse this initial condition is equivalent to the factorized one H2.17II2.41JI , We use 
the former one since it is more convenient when dealing with the general restrictions on the work-extraction. 

2. Cyclic external fields. For the following derivation, the Hamiltonian Hp(t) of external fields acting on the spin 
is completely arbitrary. In particular, it need not be composed by pulses, where it would vanish outside of the 
pulses. The only general assumption made on H F (t) is that its action is cyclic at some final time if: 

H F (0) = H F (tf) = 0. (5.1) 



We can find the following two relations (derived explicitly in Appendix iDl) 



ir_-(i-^)A^. (5.2) 

T 



W> ( 1-^ ) (Afli + Affs), (5.3) 
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where 

A# k = tr(# fc [p(tf)-p(0)]), k = S,I,B, (5.4) 

are the changes of the corresponding average energies of the spin, bath and interaction, with p(t{) being the complete 
density matrix of the spin and bath at time r, and where the total work reads: 

W = AH s + AHi + AH B . (5.5) 

Here are implications of Ecis. I|5.2I5.3[) . 

• If T$ > T and work is extracted, W < 0, l|5.2(l implies 

AH S < 0, AHi + AH B > : (5.6) 

the system looses energy, while the bath gains it and the amount of the extracted work \W\ is then bounded 
from above by |Affs|. 

• If T = Ts, both Eas. l|5.2|l5.3[l produce: W > 0, which is, in fact, the statement of the second law in Thomson's 
formulation: no work can be extracted from an equilibrium system by means of cyclic perturbations. 

• If Ts < T, inequalities in Ea. (|5.6(l are reversed: now work-extraction implies that 

AH S > 0, AHj + AH B < 0, (5.7) 

and \W\ is then bounded from above by \AHi + AHb\. 

These conclusions are close to what one could have expected from the standard (phcnomcnological) thermodynam- 
ical reasoning Q. However, it should be emphasized that in contrast to typical textbook derivations, Eqs. Q5.2I 15.3(1 
were derived starting from first principles (see Appendix[DJ), and, moreover, their derivation is by no means restricted 
to a weak bath-spin coupling, a condition which need not be satisfied in practice. 



Efficiency and Heat. 



Another useful notion is the efficiency r/ of the work-extraction, which shows how economically non-equilibrium, 
two-temperature resource is employed in work-extraction 0,0,13- The special importance of efficiency is related to the 
fact that in the standard thermodynamics it is bounded from above by Carnot's value, which is a system-independent 
quantity. 

Though our system starts out of equilibrium due to different initial temperatures of the spin and the bath, the 
notion of efficiency should be studied for it anew, since it does not automatically fall into the class of heat-engine 
models, as studied in textbooks of thermodynamics and statistical physics |!l |2j, |3j : 

• There is no working body which operates cyclically between two thermal baths. With us cyclic processes are 
defined with respect to the work source. 



The interaction between the systems having different temperatures 
the bath — need not be weak. 



-in the case discussed here, the spin and 



• We do not require that our systems always stay very close to equilibrium. In contrast, both during and im- 
mediately after the work-extraction process, the spin is in a non-equilibrium state, which in general cannot be 
described in terms of a time-dependent temperature. 

However, in spite of all these differences we can define the notion of effciency and this will be an equally useful 
characterization of the work-extraction process 0, 0, 13 ■ 

Recall that external fields are acting exclusively on the spin variables and not on those of the bath. This implies 
that when during work-extraction the source of work receives energy \ W\, this energy consists of a contribution coming 
directly from the spin and of a part which comes to the work-source from the bath but through the spin. In this 
context one can write the change of energy of the spin as 



dt 



tr(psW^sW) 



tr 



Tt Ps(t) ) H s (t) 



ti- 



ps (t) [ g- t H S (t) 



(5.8) 
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where in our case the Hamiltonian of the spin reads from Eqs. Q2.2I 13 . 1|> (note analogy with l|3.2(l 'l: 

Hs(t) = ^z+\ E h k(t)*k- (5.9) 

k— x,y : z 

The partial time-derivative in l|5.8|) stresses that we are in Schrodinger representation. When deriving (|5.8|) we have 
used dtHs(t) — dtHp(t) and H3.14|l . The last equality in (|5.8|) serves as a definition of heat (dQ) 8 . 
Integrating this from to r and using l|5.1|l and l|5.5|l we obtain 

AQ = - (AHj + AH B ) . (5.10) 

Note that in the above definition of heat, the average interaction energy is attributed to the heat received from 
the bath altough it by itself depends also on the variables of the spin; see Eq. (|2.5|) . The reason for this asymetry is 
clearly contained in the very initial statement of the problem, where we — quite in accordance with the usual practice 
of statistical physics — restricted the work source to act only on the spin. 

All this being said, one can now proceed for W < (work-extraction) with the usual definition of efficiency as the 
ratio of the useful energy \ W\ to the maximal energy involved in the work-extraction: 

„= m f511) 

/_ max(|Aff s |, \AH l + AH B \)- 1 ' 1 

For T S >T Eqs. l|B31l5TH|l and W < imply \W\ = \AH S \ - \AHj + AH B \, and then l|5TT|l results in 

-a- 



Analogously, for Ts < T we have 

\W\ _ \W\ 



(5.13) 



' \AH l + AH B \ |V^| + |AFs| 

from \W\ = \AHi + AH B \ - \AH S \. 

It is now seen from Eqs. (|5.21 15.31 15.61 15.7fl that the efficiency is always bounded by the Carnot value: 

nun(T,r s ) 
max ( T, Ts ) 



6. WORK-EXTRACTION VIA TWO PULSES. 



1. Setup of pulsing. 

Let us now dctailzc the setup of work-extraction. The spin and the bath are prepared in the state i|2.40fl with different 
temperatures Ts and T for the spin and the bath, respectively. Thus, the initial average population difference (a z ) is 
given by (|2~4H) . 

Alternatively, we can prepare the spin+bath in the state l|2.17llT4T|l . In this case, one waits for a time t ^> 1/T for 
ensuring the robustness of the results. Then the setup does not depend on details of the initial preparation, because 
the initial conditions (|2.17l IT4T|) and (|2.4U|I have become equivalent. 

The final ingredient of the setup are pulses 

V t =Vl, V t +r=V 2l (6.1) 

applied at times t and t + T, respectively. 



Note that in the equilibrium thermodynamics people frequently distinguish functions of a quasi-equilibrium thermodynamical process 
(the one which can be viewed as a chain of equilibrium states) from functions of the state. In this context the change in heat is written 
as dQ. Here we consider (possibly strongly) non-equilibrium situations, where almost any quantity (e.g. energy) is a function of the 
process. Therefore, we do not introduce the symbol d. 
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Formulas for work. 



The work done for the first pulse reads (as defined in Eas. l|3.16|) ) 

£ . „ „ . h 

- {V 1 a z - a z ) t + - 



Wx = - (Vxa z - a z ) t + - ( (Vtf z ~ a z ) X ) , (6.2) 



where for any operator A the average 

[A) =tr[ A p(t)} (6.3) 



refers to the time t just before the application of the pulse. The value of W\ is worked out by recalling the parametriza- 
tion 13.20fl . the evolution of the collective bath coordinate X as given by Ea. l|2.9|l . and finally the initial condition 
(|2~T71l2~4Hl . The final result reads: 



Wi = (!-<&>) 



(6.4) 



with G as defined in Eq. I|C5|) . and where c^l is the corresponding parametrization coefficient of the first pulse as 
defined by (Oty . 

As follows from Ts > and (a z ) < (see (I2.42f) '). the work W\ is always positive. This is in agreement with the 
thermodynamical wisdom of local equilibrium: the second term in the RHS of Eq. I|6.4|l is the contribution from the 
spin energy and it is positive, since the spin was in equilibrium before the application of the first pulse. Another 
positive term |(1 — c z J z )HG in the RHS of Eq. (|6.4|l comes from the interaction Hamiltonian (the bath operators, 
and thus the bath Hamiltonian, are not influenced by this first pulse). Again, it is intuitively expected that the 
interaction Hamiltonian should make the average energy costs higher. 

The work done for the second pulse reads analogously to Eq. Ij6.2|l : 

£ IT, " - ^ 

~ \P 2 a z - a z ) t+T + - v . ^ . 

Z Z \ / t+T 



W 2 = -(V 2 a z -a z ) t+T + -((V 2 a z -a z )x)^ , (6.5) 



where the averages (...)t+ T refer to the time just before the application of the second pulse. 
Eq. (|6.5|) is worked out in Appendix IU1 with the result for the total work W = W\ + W 2 being: 

W= -\ (l-4 2 )c«) (fr,)+«-«W»{c?i c[%e^ <e*»»-*,>} (6.6) 
'' ° (1 - c«) + \ (1 - 4 2 1) (G(r) + 92 (r) c«) (6.7) 



2 

-e 



c(T) »{cfi fle* (»R^(r) (e^^a z >+S 52 (r) (e te ^>)}. (6.8) 

The detailed explanation of various terms in this expression and of their physical meaning comes as follows. 
The first term in the RHS of Eq. I|6.6|l is the contribution from the initial spin energy. The second term comes from 
the transversal degrees of freedom excited by the first pulse. The factor e~^ T " ) accounts for the reduction of these 

(1 2) (1 2) 

terms in the time interval r between pulses. Recall that the parametrization coefficients c + ' z and c z ,' z for the first 
and the second pulse are defined in Eq. I|3.19|l : 

The terms in Eqs. (|6.71 16. 8|) are the joint contribution from the bath Hamiltonian (|2.3I) and from the interaction 
Hamiltonian l|2.5|l . The last of them couples to the transversal degrees of the spin, as reflected by the presence of 
e~^ T \ Recall that the averages (...) in Eqs. I|6.61 16. 8|) refer to the initial state Eqs. H2.17lE£lT)l , Finally, the factors 



g 2 (r)=G-G(r) = - / (6.9) 



2L 

i + r 2 r 2 

XaCr) = -7arctan(rr), (6.10) 

(the lower index 2 refers to the two-pulse situation) come from the backreaction of the spin to the bath. 

Next we note that the behavior of W = W\ + W 2 is controlled by five dimensionless parameters (see Appendix |UJ), 
which for the ohmic case reads 

W = W 1 + W 2 = ^fw(^, 7 ,^, {&,), tt\ . (6.11) 



Note that the spin temperature Tg enters only through the initial (at t = 0) (<r z ) as given by Eq. 

There are two situations within the present setup, where work-extraction is not prohibited: T > T$ and T < Tg. We 
deal with them separately, since for these cases the work-extraction effect exists in different ranges of the parameters. 
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3. Work extraction for T > T s . 



1 . 5 

1 

. 5 

w 



-0 . 5 
-1 

0.2 0.4 0.6 0.8 1 

ZT 

FIG. 6.1: Dimensionless total work w (see Eq. 16.111 in the text) versus dimensionless time tT (the waiting time between 
the two pulses) in the regime T > Ts- We compare the extracted work for different values of the initial polarization (or 
equivalently, of the initial temperature) of the spin . = 10, 7 = 1, A = 0.01, (tr z ) = —0.8, —0.5, —0.4, —0.3 (from bottom to 
top). The two pulses are given byEqs. irTHl l6~To1 KTZOl . Work-extraction disappears for larger (& z ), that is, for closer (initial) 
temperatures of the spin and the bath. 
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FIG. 6.2: The ratio |£ = = (see Eq. PTl in the text) versus the dimensionless time rr for two pulses in the regime T > Ts . 
We compare the extracted work for different values of the dimensionless bath-spin coupling constant 7. = 10, = 0.01, 
(o z ) = —0.8 and 7 = 4 (upper solid curve), 7 = 2 (lower solid curve), 7 = 0.5 (bold curve), 7 = 0.1 (dotted curve), The two 
pulses are given by Eqs. I|6.12l 16.161 10.201 . It is seen that the maximal extracted work is a non-monotonous function of the 
dimensionless coupling constant 7. 

It was seen above that the first pulse always costs work, since it is applied on the spin whose state is (initially) in 
local equilibrium at temperature Ts. However, the first pulse can do more than simply wasting work. Consider, for 
example, a ^ pulse in the y-direction 9 : 



Vl = v {r y )' (6 ' 12) 

V(ip;y) a z = e iv *» /2 a z e"^/ 2 = a z cos ip - a x sin 99, (6.13) 

V{(p\v) *x = e^' 2 axe-^v' 2 = a z sin<p + a x cos<p, (6.14) 

V iwy) °y = °y (6.15) 



This pulses excites the transversal component (a x ) which starts to decay under action of the bath, and thus correlations 
between the spin and the bath are established. The proper second pulse is then applied at time r, for instance +? in 





Here and after, we do not present the results of the full optimization of the work over the parameters of the involved pulses. The reason 
is that we do not want to make the pulsing setup too much dependent on the details of the model. On the other hand, the results of 
this full optimization did not show any qualitative difference with the presented ones. 
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the 2-direction: 



where 



V 2 = V 



(|;x), (6.16) 



V(cp;x) & z = e i<p**/2fr xe -i<p&*/2 = a z cos (p + ay sin (p, (6.17) 

V(ip;x) fry = e iv&*/2d- ye -i<p»*/2 = -fr z shxip + fr y cosip, (6.18) 

V (tp; x) fr x = fr x . (6.19) 
Note that our choice of pulses corresponds to 

<£i = l, = c«=0, 4 2 )=0. (6.20) 

It appears that not only some work is extracted by the second pulse, but the overall work by the two pulses can be 
negative for properly chosen time r: 

W = Wi + W 2 < 0, (6.21) 

as seen in Figs. 16.11 16.21 This is one of the central results of this paper. 

The time r needed for work-extraction should be neither too short — otherwise the two pulses will effectively sum 
into one, and we know that no work-extraction is achieved by a single pulse, — nor too long, otherwise the transveral 
degree of freedom excited by the first pulse will decay, and we will have two isolated single pulses. This is seen in 
Figs. 16.11 16.21 Note that the choice of pulses is obviously important for having work-extraction. Eqs. ((6.121 16.16(1 
represent only one particular example leading to work-extraction in the regime T > Tg . 

As for the magnitude of the extracted work, one notes from Eq. 16.11fl and Fig. 16.11 that it is of order of HT/2, 
which is basically the response energy of the bath. This is not occasional, since as seen from (|5.7|l . the work in this 
regime T > Tg is coming from the bath. 

Noting the ratio e/(HT) — 0.01 in Fig. 16.11 — this and even smaller ratios are usual for the realizations of the model 
as we discussed in section 0] — we conclude that the extracted work can be of several orders of magnitude larger 
than the energy of the spin. On the other hand, the extracted work is limited by ~ T which is the characteristic 
thermal energy available in the bath. Indeed, as seen in Fig. 16.11 the etracted work can be of order of HT, while the 
bath temperature is nearly ten times larger: T = lOfiF. Not unexpectedly, work-extraction disappears when the 
temperatures T and Tg are close to each other; see Fig. 16.11 

Let us return once again to the optimal time- interval r. As Figs. l6.lll6.2l show. the value of r at which the extracted 
work is maximal is roughly of the same order of magnitude as 1/T. However, the optimal r can be much larger (e.g., 
~ 10 3 /r) for smaller coupling constants 7, that is, one can increase the waiting time between the pulses at the expense 
of reducing the magnitude oc of the extracted work. 



4. Work extraction for T < Ts. 



Let us now turn to scenarios of work-extraction in the regime Tg > T. As seen from (|5.6|) . if there is work-extraction 
at all in this regime, the work should come from the average energy difference of the spin, while AHi + AHb 
is then necessarily positive. Since the latter quantity is of order of 7r (response energy of the bath), and the 
spin's energy difference is obviously of order e, there are two ways to try to achieve work-extraction, that is, to get 
W = \AHi + AHb\ — |AiJs| < 0: One should either take s/(HT) ~ 1 or take the dimensionless coupling constant 7 very 
small. The second way did not lead to work-extraction, since the required coupling constants are so small that the 
spin effectively decouples from the bath. In contrast, the first case with e/(hT) ~ 1 led to a sizable work-extraction, 
as see in Fig. 16.31 Recall in this context that systems with e/(HT) = Q/T ~ 1 are well-known; see section0]for details. 

As compared to the previous regime, here the choice of pulses has to be different for the work-extraction to be 
possible. For example, 

p 1 = v(-~;x), V2=v(~;y), (6.22) 

for the first and the second pulses respectively; see Eqs. 1(6.151 l6~1H)) for the definitions of pulses. We see from ((3.20(1 
that this choice amounts to substituting 

c +,z = *j c i,+ = 2' = C L1 — 0j (6.23) 

into Eqs. lf5iniOl6~gJl. 
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FIG. 6.3: Dimensionless work w (see 16.111 ) versus dimensionless time tT for two pulses in the regime T < Ts in the case of 
7 = 0.1, (& x ) = -0.01. The choice of pulses is given by Eqs. KTI\ l6~^3l . Full line: ^ = 0.1, ^ = 3, Dashed line: ^ = 0.1, 
# = 2- Bold line: X = 1, = 3. 



5. Efficiency of work extraction. 
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FIG. 6.4: Efficiency r] versus dimensionless time rF for two pulses in the regime T > Ts. = 10, 7 = 1, ^ = 0.01, 
(ov) = —0.8. The two pulses are given by Eqs. 16.121 16.161 16. 2011 . The efficiency is slightly below than the corresponding 
Carnot's value and is maximized over rr almost simultaneously with the dimensionless work w; see Fig. 16.11 
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FIG. 6.5: Efficiency r\ (upper curve) and dimensionless work w (lower curve) versus dimensionless time rF for two pulses in the 
regime T < T s . £ = 0.1, 7 = 0.1, ^ = 3, (a z ) = -0.01. The two pulses are given by Eqs. I|6.22l 16.23(1 . The efficiency is below 
than the corresponding Carnot's value 0.99 and is maximized over rr almost simultaneously with the dimensionless work w. 

We shall now discuss the efficiency of work-extraction as defined by Eqs. (|5.11l 15.121 15.13|) . To calculate it one 
needs to know the total work given by Eqs. (|6.4II6.6H6.8|I . and the contribution AHg to the work W coming from the 
average energy of the spin, which is read from the RHS of Eq. 1)6. 6|) . 

The efficiency as a function of rr is presented by Figs. 16.41 16.51 for T > T$ and T < Ts, respectively. There are 
several important things to note. 

• For T > Tg the efficiency can be very close to unity, if the temperatures T and Tg are sufficiently separated from 
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each other, which is the case in Fig. 16.41 It is, however, always limited by Carnot's value, as given by Eq. I|5.14|l . 
For T < Tg the effiency is sizable, but is rather below the corresponding Carnot value. 

• The work and efficiency are maximized over tT simultaneously. 

• Recall in this context that in the standard thermodynamics efficiencies close to the optimal value are connected 
to very small work per unit of time (zero power of work), since they are achievable for very slow processes. This 
is not the case with the presented setup. As seen from Figs. 16. ll 16.21 l6~31 the work is extracted on times which 
are of order of 1/r (response time of the bath), which is typicaly much smaller than the internal characteristic 
time of the spin. Thus, in Fig. 16. 41 we have nearly optimal efficiencies together with the maximal work and 
a finite power of work. 



7. WORK-EXTRACTION VIA SPIN-ECHO PULSES. 



So far we assumed that we deal either with a single spin coupled to the bath, or, equivalently, with an ensemble of 
identical non- interacting spins each coupled with its own bath 10 . However, many experiments — especially in NMR 
physics — are done on ensembles of non-interacting spins which are not in identical environment. The difference lies 
in the different energies e. This can be caused by inhomogeneous fields contributing into energy e, or by action of 
environment, e.g., chemical shifts for nuclear spins [8ft l37l l3Sj| or effective g-factors for electronic spins in a quantum 
dot. It is customary to regard these energies as random quantities, so that the collective outcomes from such ensembles 
are obtained by averaging over e = Ml the corresponding expressions for a single spin. We shall assume that the 
distribution of f2 is gaussian with average f2o and dispersion d: 11 

P(fi) = 1 e -(^-o ) 2 /(2rf). (7.!) 



It is now clear that the averaging over P(Q) the oscillating terms e lSlT will produce ~ e dT / 2 that is, a strong 
decay on characteristic times 

T* oc 1/Vd. (7.2) 

For t/T£ ^> 1 all the terms containing e will be zero after averaging, and the corresponding averaged work for 
two pulses will always be positive as seen from Eqs. (|6.41 16.61 16.71 16. 8| . Indeed, all possible negative values of the full 
work W were related to transversal degrees of freedom excited by the first pulse. These terms come with the factor 
e lflT which is connected to the free evolution in the time-interval r between the two pulses. Due to the decay of these 
terms after 3> 1, it is impossible to extract work from this ensemble via two pulses. 

However, we can extract work even in the strongly-disordered situation with T£ being short, if we combine our 
work-extraction setup with the spin-echo phenomenon pjil l35| . For our present purposes this amounts to applying a 
7r-pulse, for instance in ^-direction: 

V* o z — — <7 2 , V 7T o-y — ~d-y 1 Vtt6- x ~6- Xi (7.3) 

right in the middle of two pulses V\ and Vi (to be tuned later on) applied at times t and t + 2t, respectively. The 
work done by the first pulse reads from Eq. (|6.4(l after averaging over P(f2) given by Eq. i|7.1[l : 



(7.4) 



Wi = (l-c#) 
where 

E = --J dfiP(Q)fitanh^— < 0, (7.5) 



The assumption that each spin has its bath is a natural one for cases when the spins are sufficiently well separated \ '>(l lilt . This 
assumption is in a sense also a pess imistic one, since admitting a single bath for all the involved spins — a situation which has its own 
relevance in NMR/ESR physics l3Sl — we may get additional, collective channels of work-extraction. 

The assumption on the gaussian character of this distribution can be motivated by the central limit theorem, where the randomness of 
f! is viewed to be caused by many (nearly) independent small random factors. 
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is the average initial energy of the ensemble of spins. The work done by the 7r-pulse at time t + r is found from 16.5(1 

(2) (2) 

by substituting there the parameters c z .' z = —1 and c z ' + = of this pulse: 

W T = hG(r) + hg 2 (T) c« - 2£c«, (7.6) 

where g 2 (r) is defined in Eqs. 1)6. 9|l . It is seen that W v > 0, because the 7r-pulse does not couple properly with the 
transversal degrees of freedom excited by the furst pulse. Thus, both pulses V\ and V n waste work. 

Ultimately, the total work W — W\ + + W 2 done by the three pulses together is derived in Appendix ITU to be 



W 



HG (l + 4 2 1 c«) + hG(r) (2 - cfl - c«) - (l + ci 2 ] c« - c« - 4 2 ]) (7.7) 

4«(r)+«(2r )3? | c W c (2| ( [2 ft£( T ) _ ft £( 2r ) ] [ sin X 3 + i m cos X 3 ] - tigs [cos X 3 - i ro sin Xs ] ) } (7.8) 
4 2 14!1) + e- 4 «M + ^5R{ C W 4 2 i(2£;cosx3" i ft«osin X3 )}, (7.9) 



2 



1 



where 

g 3 (r) = G-G(2r)= i+ 7 4 ^ r2 , (7.10) 

X3 (r) = 2F(r) - F(2r) = 7 [ arctan(2rr) - 2 arctan(-jT) ] , (7.11) 
are the backreaction factors for the considered setup of pulses, and where 

to = - / dQ P(fi) tanh < (7.12) 



is the average magnetization of the ensemble. 

As compared to Eqs. (|6. 61 16.71 16.81) which present the work for two pulses, Eqs. 1)7.71 17.81 17. 9fl are different in several 
aspects. 

1. There are no oscillating factors e which after averaging over the distribution P(£l) would produce damping 
on times T 2 * . This is due to the 7r-pulse (|7.3() in the middle of two pulses (spin-echo setup). A simple explanation 
on why the terms oc e lVtT are absent is as follows. Assume that the interaction with the bath is absent and the 
spin moves under dynamics generated by the free Hamiltonian H$ = ^a z . Denote by the corresponding 



Heisenberg evolution operator: £^ A = exp 



A exp 



It is now seen with help of (|7.3[1 that the 



factor e drops out (as if the time had been reversed) : 

£(0) <p 2 E (0) &+ = e iSlr £ {0) V2 a+ = e «*rg(0) & _ = e iQx e -iQr ~_ = & _ (? 13) 

2. The decay (decoherence) factor e -4 ^ r ) + £( 2T ) in Eqs. (|7.81 17. 9| is different from e~^ 2r \ The last decay factor 
is the one generated by the free (unpulsed) evolution during the time 2r. Only in the exponential regime 

~ t/T 2 we shall have e~ 4 ^ T " )+ ^ 2r " ) ~ e~^ 2r \ (Recall that the exponential regime is present for the ohmic 
spectrum at long times, see section EOl) For gaussian decay ~ t 2 /T 2 2 , e~^ 2r ) predicts sizable decay in 
contrast to e~ Ac ^ T ^ + ^ 2T ^ ~ 1. This partial inhibition of decay due to 7r-pulse(s) is known in NMR physics j3|| 
and has been recently reinterpreted as a quantum error corection scheme 23, 24] . 

3. Now there are two independent parameters which characterize the initial state of the ensemble of spins: E and 
to. The work W in Eqs. I|7.71 17.81 17.9fl can be expressed in the dimensionless form similar to Eq. (|6.11|l : 

w ^ ( T V T s d \ 



2 \hr' t ' hV r 2 ' 

It is now more convenient to account for the temperature of the spin via , and there is a new dimensionless 
parameter which quantifies the ratio of the response time 1/T to T 2 * = 1/ Vd. The average magnetization to 
is expressed via ^ and ^ . 

Fig. 17. II describes a scenario of work-extraction in the regime Tg > T and for pulses 

Pi = v(^;x), V 2 = v(~;y). (7.15) 
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FIG. 7.1: Dimensionless work for three spin-echo pulses, = 10, 5, 1, 0.5 (from top to bottom), S = 10 3 , = 10 2 , 7 = 0.1, 
^ = 8. Work-extraction is poor or disappears for smaller or ^P-, because there is too much random thermal energy in the 
ensemble. 

This choice of pulses amounts to substituting 

c% = i, c£ 2 i = i ci 2 )=0, C W=0, (7.16) 

in C2IEEIE3- 

Recall that in the regime T$ > T there is a positive contribution to the total work coming from the bath, and 
sizable average frequencies > 5 are needed to overcome this contribution, as seen from Fig. 17. II This restriction on 
the (average) frequency is similar to the one present in the two-pulse work-extraction scenario for the non-disordered 
ensemble of spins in the regime Ts > T. 

It is seen from Fig. 17. II that the initial high-temperature ensemble of spins is strongly disordered: j4 = 10 2 3> 1. 
This ratio cannot be much larger, since there will be too much random energy in the ensemble, that is the positive 
term — - E in Eq. (|7.9I) will be too large and cannot be compensated by potential negative terms. Simultaneously, the 
average magnetization \m\ will be too small. For the same reasons there are no interesting scenarios of work-extraction 
for strongly disordered ensemble in the regime Ts < T: the average magnetization |m| is too small. 

8. CONCLUSION. 

This paper describes several related scenarios of work-extraction based on the spin-boson model: spin-i interacting 
with external sources of work and coupled to a thermal bath of bosons. The work-sources act only on the spin, since 
the bath is viewed as something out of any direct access. The model has two basic characteristic features. First, the 
transversal relaxations time Ti is assumed to be much shorter than the longitudinal relaxation time T\ . This condition 
allows the notion of local equilibrium, because once transversal components decay at time T2, the spin can be described 
via a temperature different from the one of the bath. Second, the external fields are acting in the regime of short 
and strong pulses. This feature makes the analytical treatment feasible. Both these idealizations are well known in 
NMR/ESR physics and related fields, and were applied and discussed extensively in literature |U 0, HE HE 113 • It 
may be of interest to see in future how precisely finite Ti-times and finite pulsing-times influence the work-extraction 
effect. 

The work is extracted from an initial local-equilibrium state of the spin at temperature Ts which is not equal to the 
temperature T of the equilibrium bath. As we recalled several times, Thomson's formulation of the second l aw p rohibits 
work-extraction via cyclic processes from an equilibrium state of the entire system: T = Ts |2(| I27L l28l I29L I30L l31j . 
In this spirit one would expect that work-extraction is also absent when external fields are acting only on the spin in 
a local equilibrium state pM l3S| . We have shown however, that this is not the case. It is possible to extract work in 
this latter setup due to the common action of the following factors: i) backreaction of the spin to the thermal bath; 
ii) generation of coherences (i.e., transversal components of the spin) during the work extraction process. 

With help of the spin-echo phenomenon it is possible to extract work from a disordered ensemble of spins having 
random frequencies. This ensemble can even be strongly disordered in the sense that the relaxation time 7^* induced 
by the disorder is much smaller than both the T2-time and the response time of the bath. 

As to provide further perspectives on the obtained results, let us discuss them in two related contextes, those of 
lasing without inversion and quantum heat engines. 



21 



1. Comparing with lasing without inversion. 

As we discussed in the introduction, besides the standard lasing effect, where work is extracted from a spin having 
population inversion (i.e. having a negative temperature), there are schemes of lasing which operate with a weaker 
form of non-equilibrium, since they employ three or higher-level atoms which are initially in a state with non-zero 
coherences (i.e. non-zero off-diagonal elements of the density matrix in the energy representation). There are numerous 
works both theoretical and experimental, partially reviewed in [3i|, showing that in such systems one can have 
various scenarios of lasing without inversions in populations of atomic energy levels. In quantum optics lasers without 
inversion are expected to have several advantages over the ones with inversion. 

The effects described by us also qualify as lasing without inversion (or more precisely gain or work-extraction 
without inversion). There are, however, several important differences as compared to the known mechanisms. 

• We do not require coherences present in the initial state. Our mechanism operates starting from initial local 
equilibrium state of the spin, which by itself is stable with respect to decoherence (i.e., to both T 2 and 7^* time- 
scales). It does employ coherences however, but they are generated in the course of the work-extraction process, 
which, in particular, means that all the energy costs needed for their creation are included in the extracted 
work. 

• We do not need to have three-level systems: the effect is seen already for two-level ones. 

• In one of our scenarios the extracted work comes from the bath if its temperature is higher than the initial 
temperature of the spin. Due to this fact, the extracted work can be much larger than the energy change of 
the spin. Thus the work extracted per cycle of operation can be much larger than for the standard lasing 
mechanism, where it is of order of the spin's energy. 

• The work is extracted due an initial difference between the temperature of the spin and that of the bath. This 
difference can be created, e.g., by cooling or heating up the bath on times shorter than the 7i-time. Alternatively, 
one can cool or heat up the spin with the same restriction on the times. The latter preparation of initially non- 
equilibrium state is similar to the analogous one in the standard lasing mechanism, except that no population 
inversion has to be created (i.e., no overcooling of the spin), and the spin's temperature can be increased or 
decreased. 



2. Comparing with quantum heat engines. 

The standard thermodynamic model of a heat engine is a system (working body)operating cyclically between two 
thermal baths at different temperature and delivering work to an external source [ij, |3| . The work pruduced during 
a cycle, as well as the effciency of the production, depend on the details of the operation. The upper bound on the 
cffcicncy is given by Carnot expression, which is system- independent (universal). This efficiency is reached for the 
Carnot cycle during very slow (slower than all the characteristic relaxation times) and therefore reversible mode of 
operation Q, Though Carnot's cycle illustrates the best efficiency ever attainable, it is rather poor as a model for 
a real engine. This is explicitly caused by the very long duration of its cycle: the work produced in a unit of time is 
very small (zero power). This problem initiated the field of finite-time thermodynamics which studies, in particular, 
how precisely the efficiency is to be sacrified so as to reach a finite power of work . 

In a similar spirit a number of researchers transferred these ideas into quantum domain designing models for engines 
where the basic setup of the classic heat engine is retained, while the working body operating between the baths is 
quantum Hi, US El [l|. 

Our setup for work-extraction can also viewed as model for a quantum engine. It is, however, of a nonstandard 
type since there is no working body operating between two different-temperature systems (in our case these are the 
bosonic thermal bath and the ensemble of spins). The two systems couple directly and the work-source is acting on 
only one of them. In spite of this difference, the notion of efficiency can be defined along the standard lines, and it is 
equally useful as the standard one; in particular, it is always bound from above by the Carnot value. We have shown 
that the efficiency can approach this value at the same time as the extracted work approaches its maximum. This is a 
necessary condition for a large efficiency to be useful in practice. Moreover, the whole process of work-extraction takes 
a finite time of order of the response time of the bosonic bath, which is actually much smaller than relaxation times 
of the spin. Thus, the three desired objectives can be achieved simultaneously: maximal work, maxmal efficiency and 
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a large power of work. 
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APPENDIX A: QUANTUM NOISE GENERATED BY OHMIC BATH. 

Here we discuss properties of the function: 

/>OC 

Kit) = 7 / dojuj coth(hu)(3/2) e" w/r cosut. 



-huj(3n 



In the given integration domain one can use 

coth(JM?/2) = 1 + 2 -pjrzn = 1 + 2 E ' 

n— 1 

and get from (|A1|) : 
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With help of a standard relation: 
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where ip(z) = Y'(z)/T(z), one obtains 
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= -[<//(! + K-it) + ip'(l + k + it)]. 
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Combining (|A5|) with l|A3() and k = l/{h(3T) one ends up with the following formula 
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Let us now consider separately the cases of low and high temperatures. For UTj3 1 one uses the known relation 
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and obtains from 1A6 



K(t) = 7r 
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(l + r 2 t 2 ) 2 i 2 ft 2 sinh 2 [7rt/(ft/3)]' 
For small t (t -C 1/r) -R'(i) it is positive as it should be: 



= 7 r 2 + 
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In contrast for t ~ ft/3 3> 1/r it becomes negative, namely the noise is anticorrelated, 

1 7 T 2 TT 2 1 
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At the end it is again correlated in the limit of very large times t ^> h(3 where the first term in the r.h.s. of Eq. (|A10|) 
dominates (this domain is shrunk for low temperatures). 

In the high-temperature limit hfiT <C 1 one can use in Eq. (|A6|) the Stirling formula: 



1 1 

z + 2^ 



z > 1 



and then the quasiclassical limit for the quantum noise reads (after some more simplifications) 

1 - T 2 t 2 2 7 Tr 1 



K(t) = ~fT 2 



(i + r 2 t 2 



i + t 2 r 2 ■ 



(All) 



(A12) 



In the purely classical limit the first term in the r.h.s. can be neglected and we return (for tT ^S> 1) to the classical 
white noise with the strength . 

Finally in the context of Eq. (|A6() we notice the following useful relations: 
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which are used in the main text. 



APPENDIX B 

Here we shortly outline how the two-temperature state l|2.40|l can be prepared starting from the overall equlibrium 
state 



P(0) 



1 



■ exp 



-/3s#s -/%(-»! + #b) 



(Bl) 



which has equal temperatures of the spin and the bath. 

Assume that the bath was subjected to another much larger thermal bath (superbath) at temperature T different 
from Tg, so that the total Hamiltonian of the spin, bath and superbath reads: 



-Htotal = H + H s 



sup, 



(B2) 



where the latter operator H sup characterizes the weak interaction of the bath with the superbath and contains also 
sclf-Hamiltonian of the superbath. Thus 



[-ffs, H sup ] — 0. 



(B3) 



Now the statement of this appendix is that under the action of the superbath at temperature T, the common state 
of the spin and the bath will relax to the state Eq. ({Bill or 12.4011 with different temperatures for the spin and the 
bath. The reason is that due to Eq. (|B3|I . a z is conserved during the whole evolution generated by the superbath, so 
that a z does not relax and keeps its value given by Eq. (|B1|) 12 In contrast, the variables of the bath — including X— 
do not have such a protection, so they relax under influence of the superbath. Let us now substantiate this statement. 

Because [Hs,H] = 0, the initial equilibrium state p(0) of the spin and the bath can be represented as 



P(0)= PiiPiMtiM, 



(B4) 



3=±1 



12 In a more realistic situation, where Ti-time is kept finite, the relaxation time of the bath under action of a superbath has to be much 
smaller that T\ , so as to create the temperature difference between the spin and the bath. 
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where 



e -]Pse/2 

2cosh(/3 S £/2) : 



3 = ±1, 



(B5) 



are probabilities for the spin to be up or down respectively, \j) is the eigenstate of H$ = §<r 2 with eigenvalue j = ±1, 
and where 



z 



exp 



k ( |i • //„ 



= tre e 



i = ±i, 



(B6) 



are conditional states of the bath. 

The total initial state of the spin, bath and superbath thus reads: 

Ptotal(O) = PJJ Amp(O) <8> p j3 (0) \j)(j\, 



(B7) 



where /9 SU p(0) is the initial equilibrium state of the superbath. Note that due to weak coupling between the bath and 
superbath, their initial states can be assumed to be factorized. 

As follows from (|B3IIB4J| . the time-dependent state of the total system consisting of spin, bath and superbath can 
be presented as 



where Cljj(t) — the conditional joint state of the bath and superbath — satisfies the von Neumann equation 



(B8) 



Ml 



j j 



[-^X + Hb + H sup , Qjj], 



with the initial condition: 



%(0)=p sup (0)®fti(0). 



(B9) 



(BIO) 



Thus, fljj moves according to the Hamiltonian + Hb + H sup . It is now clear that in the weak coupling limit of 
the bath-superbath interaction the marginal conditional state tr sup f2jj(t) will — for sufficiently long times t — relax 

to Gibbs distribution at temperature T (equal to the one of the superbath) and with Hamiltonian + Hb- Thus 
the (unconditional) marginal state of the spin and the bath will indeed relax to 



p cx exp 



-p s H s - P{Hj + Hb) 



(Bll) 



APPENDIX C 



Here we explain in detail why the initial conditions (|2.17l |2~TTJ) and (|2.40|) are equivalent. 
One can write the full Hamiltonian H defined in l|2.1|l as 



k 

and diagonalize it via a unitary operator: 
U = exp 
Thus the operators 



t , 9kcr z 



2u!k 



a-k 



2t0k 



Ua k t/ f 



ak 



h = a k + 



9k^z 
2u k ' 



^k 



9kOz 

2u k 



Ua z tTt =a z 



[b k ,bj] = S k i 



(CI) 



(C2) 



(C3) 
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are distributed — over the initial state (|2.40(l — independently from the operators of spin. Moreover, as follows from 
(jCll I2.4U|) . the operators bk have on the state l|2.40|l exactly the same statistics (i.e., the same correlators) as the 
corresponding operators &k on the factorized state (|2.17|) . 

Now note that for the initial condition l|2.40[l . a z (0) and the quantum noise operator r)(t) are in general not 
independent variables, in contrast to the case of the factorized initial condition (|2.17l I2.41|) . However, for t ^ 1/T 
they do become independent: 



fi(t)=Ut) + vAG(t)-G), %(t) = ^ 5fe [6[(0)e^ t + 6 fc (0)e-^*], (C4) 

k 

where 

G = V 2k. (C5) 

k 

is the limit of G(t) for t ^> 1/T. Taking the latter limit in l|C4|l . one gets is equal to fjb{t) and is thus independent 
of a z . Recalling that fib{t) has on the state 1|2.4U|) the same statistics as fj(t) on the factorized state l|2.17|l . finishes 
the argument: the equivalence holds for times larger than the bath response time (1/T for the ohmic situation). 

Note that the thermodynamical limit for the bath is essential for this conclusion. Otherwise, G(t) would qbe a 
finite sum of cosines, and would not converge to G. 

APPENDIX D: DERIVATION OF EQS. J5721 15731) . 

Assume that the initial state of the spin and bath is: 

P (0) = i exp [-foHs - P(Hi + Hb)] , Z = ixe -^-^+^\ (Dl) 

with different temperatures for the spin and the bath. 
An external field V(t) is acting on the system, 

H{t) = H + V{t) (D2) 
such that it is zero both initially and at the moment t = r: 

V(t) = V(0) = 0. (D3) 
This condition defines cyclic process. The total work which was done on this system reads: 

W = AH S + AHi + AH B , (D4) 

where 

H k = tr(^H k [p(r)-p(0)]), k = S,I,B, (D5) 

are the changes of the corresponding energies, and where p(r) is the overall density matrix at time r. 
Recall that the relative entropy (see, e.g. [2(1]): 

S[p\\a] =tr (pin p- phi a) > 0, (D6) 

is non-negative for any density matrices p and a. One now uses: 

S[p(r) | \p(0)] = tr( p(r) In p(r) - p(r) In p(0) ) = tr ( p(0) In p(0) - p{r) In p(0) ) 

= (3 S AH S + /3(A#i + AH B ) > 0, (D7) 

where we used I|D1(I and ti p(r) lnp(r) = trp(0) lnp(0) is due to the unitarity of the overall dynamics generated by the 
time-dependent Hamiltonian H(t). 

Combining HD7|) with l|D4jl one gets Eas. lfOO|l 

W>(i-^-)aH s , W>(l-^]{AHj + AH B ). (D8) 
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Finally note that would we use the initial conditions 

p(0) = pb(0) ® PB (0) = <g> ^e"^ (D9) 

we would not be able to conclude from the above derivation that the efficiency is limited by the Carnot value. Indeed, 
instead of Eqs. (|D8|> one has, respectively: 

(3 S AH S + (3AH B > 0, W > AHi + (l - AH S . (D10) 

The latter inequality is not informative with respect to Carnot's bound, since it cannot and should not in general be 
excluded that AHi is sizeable. 

However, for the model studied in the present paper, the equivalence of the initial conditions I|D1|) and (|D9|I is 
known from other places. 

Let us emphasize the main points by which the present derivation differs from the standard textbook one: 

• No postulates were used: the whole derivation is based on the quantum-mechanical equations of motion and 
certain assumptions on the initial conditions. 

• It was not assumed that the interaction between the system and the bath is small, a restrictive assumptions 
which need not be satisfied in reality. 

• The fact of using the initial conditions in Eq. IjDlfl is important in the present derivation, though presumably 
Carnot's bound is valid in certain more general cases, such as, in our case, factorized initial conditions from Eq. 

APPENDIX E: SOME CORRELATION FUNCTIONS. 

In this appendix and in the following ones we study veraious dynamical aspects of the model defined by Eqs. H2.ll 
12.21 12. 312. 5|l . The initial conditions are given by (|2.17llT^T|) . Hereafter (...) means averaging over this initial condition. 
Let us define some correlation functions. 
a. For 

ta > h > h, (El) 
and recalling definitions 1(2.111 12.15(1 one derives using Wick's theorem in the same way as when deriving (|2.25() : 

fj(t 3 ) ft±(ti,ta) ) =±f^ l { 2k + l)\ I 2 -/ 2 dsi - ds2fc+1 (Tlv(t3)f)(s 1 )...f)(s 2k +i)}) 



(E2) 



= ±i dsK T (t 3 -s) exp -- / / ds 1 ds 2 K T {si - s 2 ) 
Jti L ^ Jti Jti 

= ±[ii{t s - *0 - ii(t z - h) + G(t 3 - tO - G(«3 - t 2 ) ] e -eto-ti)+iFfr-ti). (E3) 
where for deriving the last line we used the definition of Kq-(t): 

K r {t) = K{t) - iG{t) = £(t) - iF(t). (E4) 
Note that for t 3 = t 2 we can derive (|E3(1 in a simpler way by employing (|2.25|) and 

'fj(t 2 )n±(ti,t 2 )) = Tid t2 (n±(ti,t 2 )) . (E5) 



b. A correlation function (Jl±(ti,t 2 ) fj(t 3 )j under the same condition (|E1(1 is studied similarly to (|E3|) . the 

only difference being that the time-ordered correlation function Kq-{t 3 — s) in l|E2|l is substituted by the analogous 
time-antiordered one (time-antiordering comes due to {Eljl ) 

K A (t 3 - s) = K}(t 3 - s). (E6) 

These two functions are related by complex conjugation, as seen from (12.21(1 . Thus, 

(n ± (t!,t 2 ) fj(t 3 )) = ±[it(t 3 - tO - *&3 - *a) - G(t 3 - tO + G(t 3 - t 2 )] e -«*»-*0-H^(* a -*i). ( E 7) 

As compared to <|E3(I . the sign of G- factors is seen to change. 
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c. A correlation function between two Il-factors for 

U > U > t 2 > *i, 

is worked out as follows. First one notes 



TL ± (t3,t4,)IL T (t u t 2 )) = (Texp 



1 1 



±i I ds 4>(s) fj(s) 



(E8) 



(E9) 



where 



(s) =-1, h<s<t 2 , 
= 0, h < a < i 3 , 
= 1, t 3 < s < i 4 . 



(E10) 



One gets 



n±(t 3) t4)n T (ti,t 2 )) =exp 



1 

2 

r*2 r*4 



dsids 2 ^T(si - s 2 ) - - / / dsi ds 2 -K"r(si - s 2 ) 



^4 /"i4 



/>t2 /"C4 

/ / dsi ds 2 ^r(si - s 2 ] 

Jt! Jt 3 



= exp [-£ (i 2 - h) - £(t 4 - t 3 ) - S(*4 - * 2 ) + C(*4 - h) + £(i 3 - i 2 ) - £(*3 - ti)] 

x exp [*F(i a - ti) + %F(U ~ *a) + iFfa ~ *a) ~ iFfa - *x) - *-F(*3 - ^) + iF(t a - h)} . (Ell) 

APPENDIX F: EVOLUTION OF THE QUANTUM NOISE AND Il-FACTORS UNDER HEISENBERG 

DYNAMICS. 



Note from (|2.8I I2.lf)l how the quantum noise and Il-factors evolve under Heisenberg dynamics: 
St rj(r) = e ltit ' h fj(r) e"^*/* = f,{t + r) + a z [ G{r) - G(t + r)] , 



s t n± (t x , t 2 ) = e lkt ' h ri± (ti , t 2 ) e' lkt ' h = r exp 



= U ± {t + t u t + t 2 ) exp (±i u z [ F(*a) - F(ti) + F(t + t 1 )- F(t + 1 2 ) ] ) . 



(Fl) 



(F2) 



When deriving these equations we used [t)(t + 1), & z ]=0. Recall that a z is conserved under evolution generated by H: 
S t a z = a z . 



APPENDIX G: DERIVATIONS FOR TWO PULSES. 



The work done by the second pulse is defined as 
\\V 2 = \ (V 2 (H s + Hi) - {H s + Hj) 



| {V 2 a z - a z ) t+T + \ (<J> 2 a z - a z ) x) ^ (Gl) 

§ (4 2 1 - 1) (^) t+T +\{cfl-i)(a z x) t+T +n»{ c a<fi- + )t+r} + »te ( a +*) t+ Y 



(G2) 



where the averages are taken at the time t + r immediately before the second pulse, and where we used the definition 
(|3.19l) of the parametrization coefficients. For clarity we recall that definition here 



V k o a = ^2 c;j<7&' a = ±,z, k=l,2. 



(G3) 



b=±,z 
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In order to calculate W2 we thus have to determine X(t + r), <7+(i + r), and a z {t + r). Recall from Eq. I|3.12[l that, 
e-g-, 

X(t + T)=S t Vi£ T X, (G4) 
where £ t is the free evolution (super)operator defined in Eq. (|3.11|) . One infers from (|2.81 12.9(1 

£ T X = fj(r) - & z G(t), (G5) 

£ t)(t) = 77(i + r) - a z [G(r) - G(t + r)], (G6) 

and then 

X(t + t) = E t V x S T X = fj(t + t) + [G(t) - G(t + t)]& x - G(t) £ t V x o z (G7) 
The formula for a z (t + r) is more straightforward: 

a z (t + r) =£ t Via z , (G8) 

(*,(< + r))=c#<*,>, (G9) 

where we noted that in 

Pi **) - £ eg ($ frn) = cjg , fc = ±, z, (G10) 

n— ±,2 

only one term contributes, since (a±) = ( e lX(7z a±) = due to the initial conditions i|2.17ll2~4*T)l . 

In the same way one calculates 

(a z X) t+T = [G(r) - G(t + r)]c« - G(r), (Gil) 

a+(t + T ) = Et Vi £r <J+ = e- iF ^ +iClT E t Vxil+(0, t)&+ = e -«W^+>x*»n + ( t) t + T )£ t Vx &+, (G12) 
where we used (|F2|) . [<7 Z , U±] = 0, and where by definition (from 1F2H ): 

X(r,t)= I' ds[G(s)-G(t + s)} = F(t) + F(r)-F(t + T). (G13) 
Jo 

Now let us recall Eq. I2.25|l : 

( n±(ti,t 2 ) ) = exp[-£(t 2 - ti) + 1 F(t 2 - h)], (G14) 
because it is used in averaging the RHS of Eq. I|G12|) : 

(*+(* + r)> = c £> e --«M (e**- a z ), (G15) 



where we additionally employed the reasoning which led us to (|G10|) . 

The last term we have to calculate is ( a+ X ) . Directly multiplying 1G7|) and 1G12I) one gets 



t + T 

a+x) =e- tF ^ +inT+l ^fl + (t,t + r)£ t V 1 a + fi(t + T) (G16) 

/ t + T 

+ [ G(r) - G(t + t) ] e-^M-HOr-Hx*.^ ( t) 4 + r ) £ t{ £ + ) &z } ( G 17) 
-G(r) e- iF W +inT+ix **n+(i, t + r) £{7>i <7+<xJ. (G18) 

Following to JUU) we now expand "P\ <j+ in (|G17hlGT8|) . With the same reasoning as for IjGlOf) . we need to keep in 
these expansions only terms propotional to cv} z (since (a±) = according to the initial conditions i|2.17ll2~4*T|) '). After 
further simplifications with help of (|E7|I we obtain 



(a + i) t+r = 4^ T - fM (it(r) {j* a *a z ) + [G{T)-G{t + r)] (e**-)) . (G19) 
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The final formula for the work reads: 

\w 2 = 1(1- c«)(l - 4 2 1) G(r) + \ (1 - c« 4 2 1) G(t + r) + OSR {<£>. (a + ) t+r } + {eg. («r + x}^ J .(G20) 

Note that in the limit t — > oo (which means i 3> 1/r), one has G(t + r) — > G, where G is defined in (|C5|) . 

Eq. (|O20|) can be put into dimensionless form as anounced by (|6.11() . To this end note from (|A13I IaT4|) that £(r) 
and pC( T ) can be expressed via dimensionless quantities iT, 7 and T/(hT). In the same way we note from (|2.38l 
12.39(1 that jtG(t) and F(r) are expressed via 7 and rT. 

APPENDIX H: DERIVATIONS FOR THREE PULSES (SPIN-ECHO SETUP). 

Now we consider three pulses, V\, 7-V and V2 which are applied, respectively, at times t, t + r and t + 2t. The 
pulses V\ and Vi are kept arbitrary, while 7-V is the 7r-pulse defined by Eq. 1)7. 3fl . 

The work done for the pulse P2 is defined by the same formula l|G2|) . where now all the averages are taken at the 
time t + 2t immediately before the application of V2 ■ Our calculations in the following will be relatively brief, since 
in essence they follow to the pattern of calculations in the previous appendix. 

For ct + {t + 2t) we get 

a+ (t + 2r) = s t Vi £ r v« £ T a+ = e -iF(T)+i«>T St Vi Er n + ( 0) T y a _ 

= e -4^(r) +J F(2r) ^ ^ jj + ( T> 2r ) ft. (0, r) <7_ = e -4iF(r)+iF(2r) £ | jj + ^ 2t) ^ r )j ^ {-^ - _ | 

= e -UF(r)+iF(2r) f[ + ( t + T ,t + 2T)fl_(t,t + t) e^* 3 ** £ t {Vi &-} , (HI) 

where we we used ((F2() and defined 

X3 (r, t) = 2F(t) - F{2t) - 2F(t + t) + F(t) + F(t + 2r). (H2) 

Taking in this equation the limit t^$> 1/T and using (|2.39l) we return to the quantity X3( T ) as defined by l|7.11[) . 
With help of (|E11|) and the reasoning of (|G10(I one has 

(a+(t + 2r)> = c% e - 4 «(-)+«( 2 -) ( e -^ a z ) (H3) 
In the same way as for <|H1|I we have 

a z (t + 2r) = St V x £ T V, £ T a z = -£ t V x a z , (H4) 
while applying (|G5IIG6I) one derives: 

X(t + 2t) = £ t Vi £ T V„£ r X = fj(t + 2t) + [ 2G(t) - G(2t) }£ t V 1 a z + [ G(2t) - G(t + 2r) ] a z . (H5) 
The only non-trivial relation in calculating ( tr + X ) is 

\ / i+2r 

( fl+{t + r,t + 2r)fl-{t,t + T)T!{t + 2t)J) = [2^(t)-^(2t)-2G(t) + G(2t)] e - 4 «W+«( 2 ^) ; (H6) 

which is obtained in the same way as 1JE7I IE 1 1|) . The easiest way to check this relation is to follow to the derivation 
of (|E5J| . that is, to differentiate (|E11|) over £4, to put = t + 2r, £3 — t2 — t + r, t\ = t, and then to change the sign 
of all G-f actors in the final expression. 

If the reader has followed us so long, he/she can continue alone, since the remaining calculations are fairly straight- 
forward. 



